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ABSTRACT 




An experiment was performed to study the effects of an 
underwater explosion on a submerged test panel. This is a 
part of an on-going research project sponsored by the Defense 
Nuclear Agency (DNA) to help predict the dynamic responses 
of naval structural members due to underwater shocks. An 
important element of this thesis was to validate a finite 
element/finite central difference computer code developed to 
forecast shell responses. Emphasis was placed upon attaining 
stiffener tripping and collecting invaluable dynamic flat 
plate tripping data. 
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INTRODUCTION 



I . 



A. OBJECTIVES 

Current building specifications for U.S. Navy ships and 
submarines require shock to be met on shipboard equipment 
and that underwater shock testing be performed on at least 
one ship in each class [Ref. 1] . Underwater explosion 
(UNDEX) research is conducted at NPS because of the impor- 
tance of shock testing to the Navy. The objectives of this 
thesis are four-fold. First, to investigate the large de- 
flection elasto-plastic transient response characteristics of 
submerged stiffened flat plates, subjected to underwater shock. 
Secondly, to perform underwater explosion testing using a 
stiffened flat plate and to validate the predicted results 
with those of the experiment. Thirdly, to investigate stiffener 
tripping and its effect to the gross shell responses. And 
finally, to acquire the knowledge and technical expertise 
necessary to conduct these types of experiments. 

What is stiffener tripping and why is it significant to 
naval construction? Stiffener tripping is defined as the 
lateral-torsional instability of a stiffener as it becomes 
suddenly unstable and fails under load [Ref. 2] . Tripping 
is identified by a twisting of the stiffener about its line 
of attachment to a structure. The purpose of stiffeners is 
to add strength to ship hulls. Removal of the stiffener 
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weakens the structure. A report published in 1979 by John 

Adamchak stresses the importance of stiffener tripping 

effects on shipboard hulls and grillages by stating that 

previous underwater tests 

clearly demonstrate the significance of lateral- 
torsional instability (tripping) as a primary ductile 
failure mode for a ship's structure. [Ref. 2] 

Once a stiffened structure has experienced tripping, its 

ability to continue to resist deformation under loading is 

drastically reduced. 

B . BACKGROUND 

This is the third NPS thesis in the field of UNDEX and 
their dynamic effects on submerged plates. This research is 
part of a major project coordinated and sponsored by the 
Defense Nuclear Agency. During the early part of 1984 , Lt . 
Tom Rentz conducted the first NPS thesis research on UNDEX 
responses of flat plates. A plate positioned on top of an 
air-backed test platform was chosen to best simulate the 
conditions present in ship hulls. 

Fabrication and exact specifications of the test plat- 
form are covered in detail in [Ref. 3], with Figure 1.1 
showing a sketch of this structure. Flat plates are manu- 
factured out of 6061-T6 aluminum because it is accessible, 
easy to machine and this material exhibits strain rate inde- 
pendent properties. The first plate tested is shown in 
Figure 1.1 and consisted of two parallel rectangular bar 
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Figure 1.1 Test Structure Drawing 
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iy (typ) 



stiffeners with a thickness of 0.1875 inches and height of 
1.0 inches running transversely on the plate. 

In February 1984, an explosive shot was performed on the 
first stiffened plate. Last-minute changes made it necessary 
for sixteen half-pound rectangular blocks of TNT to be sub- 
stituted for a single eight-pound charge. Test results 
revealed large deformations and tears at stiffener ends 
which prevented stiffener tripping. 

Lt. Nelson King was the next to pursue the area of UNDEX 
at NPS . During the summer months of 1984, he performed two 
explosions, using different plate designs for each. His 
first plate was similar to Figure 1.1 except that the 
stiffeners were reconfigured to a tee, see Figure 1.2. It 
was believed by exposing more surface area to the shock 




Figure 1.2 



King's Tee-shaped Stiffener 






wave , a larger compression zone would be developed, encourag- 
ing tripping. A single nine and a half-pound test charge 
was exploded. Test results were similar to the first explo- 
sion. A move towards a simpler test panel was taken. 

The test panel was redesigned to a single bar stiffener 
with increased thickness to 0.25 inches and height to 1.25 
inches running longitudinally. Figure 1.3 shows the geometry 
of this new panel. A half-pound charge of TNT was exploded 
over the test plate. Recorded sensor results were invalid 
for this experiment but micrometer measurements taken later 
indicated that perhaps a slight amount of stiffener bowing 
was taking place and further testing with this test geometry 
was warranted [Ref. 4] . 

The test panel investigated for this thesis is identical 
to the single longitudinal stiffener model. A charge weight 
of eight pounds TNT, a depth of four feet and a plate stand- 
off distance of nine feet was selected. The experiment was 
conducted at West Coast Shock Facility (WCSF) at Hunter's 
Point Shipyard, California. 
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VIEW A-A 




Figure 1.3 King's Plate and Current Model 



II . 



UNDERWATER EXPLOSION THEORY 



Before continuing with the experimental aspect of this 
research, an overview of UNDEX theory is provided. Explosions 
are a very complex field of study. The destructive effect 
of an explosion is dependent upon the composition and geometry 
of the explosive element, the media or surrounding environ- 
ment, the stand-off distance to the target and the physical 
nature of the target. Simply stated, an explosion is a chemi- 
cal reaction which releases energy; both thermal and mechanical 
[ Re f . 5 ] . 




Figure 2.1 Shock Wave Pressure Profile 
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A. THE SHOCK WAVE 



An immediate result of underwater detonation is a high 
gas pressure build-up. The gases interact with the surround- 
ing water causing compression and a radial outward flow. 

This forms the shock wave which comprises about half of the 
explosion’s energy. It is the effects of the shock wave on 
the flat plate that is the primary area of interest for this 
thesis research. The exact character of the shock wave near 
the charge is not well known, thus, focus will be on the 
characteristics of the shock wave at some distance away. 
Experience has shown that the shock wave travels at the speed 
of sound in water (approximately 5000 feet per second) [Ref. 
6]. Formulas have been developed for distances between 10 
and 100 time the charge's radius to describe such features 
as: peak pressure, pressure as a function of time, and time 

decay constant. Ideally, the shock wave, which delivers a 
pressure incident, has a profile similar to Figure 2.1. There 
is a sharp peak pressure spike, followed by an exponential 
decay, with about one third the original pressure occurring 
after only a single time decay constant. 

A few equations useful in the description of the shock 
wave pressure profile are listed below. Pressure is des- 
cribed by equation 2.1. 

( t- t x ) 

P(t) = p e ; given t > t-, psi (2.1) 

max 1 
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0 is the explosive decay constant and P 

r - 1 max 

given by equations 2.2 and 2.3. 



is peak pressure 



0 




msec 



( 2 . 2 ) 



P 



max 




psi 



(2.3) 



The time of shock wave arrival at the target, T, is calcu- 
lated from equation 2.4. 



W is explosive weight in pounds 

R is stand-off in feet from charge to target 

t-t^ is elapsed time after shock's arrival 

A^, A 2 , and K 2 are constants which vary depending 
on the type of explosive (see Table I) 

C is the speed of sound in water 

The experimental values expected for an eight-pound charge of 
TNT with a nine foot target stand-off are peak pressure of 
3814 psi, decay constant of 0.153 msec and shock arrival time 
of 1.8 msec. 



T 



R 

C 



sec 



(2.4) 



where : 
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TABLE I 





Explosive 


Constants 




Explosive 


TNT 


HBX-1 


Pentolite 


K 1 


22,505 


22,347.6 


24,589 


A 

1 


1.18 


1.144 


1.194 


K 2 


.058 


.056 


.052 


A 2 


-.185 


-.247 


- .257 


K 3 


12.67 


14.14 


12.88 


K 4 


4.268 


4.761 


4.339 



In order to obtain clean data the secondary effects of 
UNDEX had to be eliminated or reduced. UNDEX secondary 
effects to a target include gas bubble, bulk cavitation, and 
surface and bottom reflection effects . A description of 
each of these phenomena is provided. 

B. THE GAS BUBBLE 

The gas bubble is an UNDEX phenomenon which is closely 
related to the shock wave. At detonation, while the shock 
wave flows radially outward, the highly compressed gases 
form a bubble which expands in size with an accompanying 
decrease of bubble gas pressure. The internal gas pressure 
gradually decays as the bubble expands until hydrostatic 
pressure finally takes over forcing the bubble to contract. 
Bubble contraction is slow at first but progresses rapidly 
until a minimum bubble radius is attained. The bubble even- 
tually collapses producing a pressure pulse which stimulates 
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radial flow outward forming a new bubble that initiates the 
cycle all over again. Pressure pulses continue to be pro- 
duced until all of the bubble energy is emitted. These 
pulses make up part of the reloading nature of an underwater 
explosion. Although the bubble pressure pulses are signifi- 
cantly weaker in magnitude than the incident shock wave (on 
the order of 10-15 percent of the first bubble pulse) , they 
still represent important dynamic loads to a target [Ref. 7] . 

An interesting aspect of gas bubble theory is the migra- 
tion of the bubble to the surface. The migration is upward 
due to the influence of gravity, with the most rapid changes 
noticed when the bubble's volume is at a minimum. Figure 2.2 
depicts the relationship between bubble expansion, pressure 
pulse emissions and bubble migration. This figure is adapted 
from [Ref. 7] . Experiments have shown that the energy asso- 
ciated with the gas bubble is about 41 percent of the total 
energy available during an explosion [Ref. 1]. The following 
empirical formulas are presented to give some appreciation 
for the gas bubble and its behavior. The maximum bubble 
radius and the period of the first pressure pulse are given 
by equations 2.5 and 2.6. 



A 



max 




ft. 



(2.5) 



T 




sec 



( 2 . 6 ) 
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Figure 2.2 Bubble Historical Representation 
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where :: 

A is maximum bubble radius 

, .. max 

L i 

K, and K. are constants from Table I 
3 4 

D is the submerged charge depth in feet 

Both of these equations employ the fact that the pressure at 
the surface of the water is one atmosphere or 33 feet of sea 
water . 

C. BULK CAVITATION ZONE 

Bulk cavitation is an area where air bubbles or pockets 
have formed which are a result of the reflection of the 
compressive shock wave off the water's surface. As the 
shock wave reflects from a free surface it changes in sign 
to a tension wave. The water mass eventually pulls apart 
beneath the free surface to prevent any further tension and 
voids are created [Ref. 6] . The parameters which affect the 
area of the bulk cavitation zone are charge type, weight and 
depth . 

D. FREE SURFACE AND BOTTOM REFLECTIONS 

Free surface and bottom reflections are termed 'late- 
time' effects. The shock wave moves radially outward from 
the charge, very much like the resultant water motion from a 
pebble thrown into a still pond. The incident shock wave 
proceeds directly to the target, while the bottom and free 
surface serve as boundaries or obstructions. Surface cut-off 
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occurs at the target when the reflected free surface negative.* 
pressure wave strikes. To help visualize what is happening 
with surface reflection, an image source is used. Figure 2.3 
shows the assumed image source, charge and point of interest/ 
target. Figure 2.4 depicts the type of pressure profile 
experienced by the target. Using the geometry of Figure 2.3, 
surface cut-off time and pressure can be computed by equations 
2.7 and 2.8. Both of these figures are taken from [Ref. 1]. 




Figure 2 . 3 Image Source 
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Figure 2.4 Incident and Surface Reflected Pressure Waves 

R. - R 

t c = — ^ — sec (2.7) 

w l/3 A 1 

p c = _K 1 ( ^ L r— > P si <2 - 8) 



where : 



R^ and R are expressed in feet 
C is speed of sound--500 fps 

The influences of the sea bed or bottom are not as well 
understood as those of the surface. However, the bottom is 
considered to be a perfectly rigid boundary reflecting waves 
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back with the same sign and magnitude. This is really an 
approximation as pointed out in [Ref. 6]. Figure 2.5 obtained 
from [Ref. 9] is used to graphically incorporate the effects 
of the shock wave, surface and bottom reflection pressure 
waves striking the point of interest. 




Figure 2.5 Shock Wave Interactions 
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E. EXPERIMENTAL CONDITIONS AND PRE-SHOT CALCULATIONS 

As mentioned earlier one of the objectives of this thesis 
was to obtain clean data, looking primarily at the effects 
of the incident shock wave. By careful attention to test 
platform positioning in the water, most of the secondary 
effects could be eliminated. A charge depth of four feet was 
selected which allowed for adequate venting of the gas 
bubble to the atmosphere prior to hitting the test plate. 
Additionally, experience has shown that a satisfactory vent- 
ing depth is one which is between 50 and 80 percent of maximum 
gas bubble radius [Ref. 1] . Using equation 2.5, a maximum 
bubble radius of 7.6 feet was computed. This represents a 
depth to radius ratio of 52.6%. 

To eliminate bulk cavitation, the cavitation region had 
to be first determined. The equations necessary to calculate 
this region are extensive and tedious, thus a computer program 
was used to graph the cavitated region, see Figure 2.6. To 
avoid bulk cavitation and in keeping with past NPS UNDEX 
experiments, a nine feet stand-off distance was used. 

Bottom and surface reflections are late time effects. 

By positioning the test structure directly beneath the charge, 
there is a time window of approximately 2.5 msec where 
reflections can be ignored. The experiments' sensors were 
sampled during this time frame. The depth of the water at 
the test site is approximately thirty feet. 
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Figure 2.6 Bulk Cavitation Zone 
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III. THE COMPUTER CODE 



The computer code selected for use in this thesis is 
EPSA, Elasto-Plastic Shell Analysis, developed by Weidlinger 
Associates. This is the second attempt at NPS to model a 
flat plate based on the EPSA code. Lt . Rentz was the first 
to utilize EPSA, see [Ref. 3] . 

EPSA is a finite element, finite central difference com- 
puter code. For calculation purposes EPSA assumes a symmetric 
loading pattern and time zero commences when the spherical 
shock wave first encounters the test plate. DNA/NAVSEA funded 
its development for the purpose of analyzing dynamic loading 
effects on submerged stiffened shells. The theory incorporated 
in EPSA is explained in the EPSA User's Manual [Ref. 8] . 

The purpose of this chapter is to describe how EPSA was 
used in the pre-shot calculations for this experiment. In 
Chapter V, EPSA will be reviewed emphasizing its post-experiment 
capabilities. EPSA is operated on a Vax-11/780 computer 
system located in Spanagel Hall. 

Symmetry conditions were incorporated to model the aluminum 
test plate for EPSA, see Figure 3.1. Requiring only half 
of the plate to be discretized saved on computer time and 
storage space. EPSA computes quantities from nodal points 
and elements. Plate displacements and pressures are nodal 
point values, whereas stresses and strains are elemental. 
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Selecting a discretization of fifteen columns by twenty five 
rows yielded rectangular block sizes of 0.6 by 0.48 inches. 

The total number of elements and nodes are 375 and 416, 
respectively. Sides one, three, and four have fixed boundary 
conditions and side two represents a symmetry edge. 

A computer data deck is necessary to run the EPSA code. 

A sample data deck created for this thesis is reproduced in 
Figure 3.2. The User's Manual describes in detail the infor- 
mation included in this input deck. It is through the input 
deck that the plate physical characteristics and the under- 
water experimental conditions are passed to the computer code. 

To facilitate in the viewing of EPSA output files, Francois 
Daube , a thesis student, wrote an interface program with 
PATRAN-G. PATRAN provides a color graphics capability for 
the display of von Mises ' stresses and plate displacements 
at requested time increments. Daube ' s work is presented in 
[Ref. 10] and the EPSA-PATRAN interface is covered in [Ref. 3]. 

A color chronological history of von Mises' stresses and 
displacements for an eight-pound TNT charge was created. 

By looking at PATRAN and EPSA, pre-shot outputs strain gage 
locations were selected. See Figure 3.3. Emphasis was 
placed upon positioning the gages in an X,Y orientation and 
maintaining symmetry in order to compare experimental load- 
ing patterns to EPSA predictions. 
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Figure 3.2 A Sample EPSA Computer Deck 
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Figure 3.3 Strain Gage Locations 
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IV. UNDERWATER SHOCK TESTING 



A. INSTRUMENTATION 

The UNDEX experiment was conducted at the West Coast 
Shock Facility (WCSF) . This is a Navy sponsored activity 
accessible to NPS , located about 120 miles to the north, 
in the South Bay of San Francisco. 

Thirteen strain gages, two accelerometers and two pressure 
transducers were installed on the test structure. The loca- 
tion of the strain gages has already been shown on Figure 3.3. 
The accelerometers were mounted alongside each other at the 
bottom of the test structure. One of the pressure transducers 
was placed on a nine foot arc from the charge approximately 
five feet away from the test platform to record incident 
pressure data. This was secured in place by light manila 
guidelines and anchored vertically by a weight. The other 
was situated on top of the plate coming out of one of the 
corners at a 45° angle towards the plate's center. This 
transducer served as the basis for the plate's total pressure 
measurements. The sensor specifications are listed in Table 
II . 

Two Honeywell M101 Wideband II (direct record) tape units 
operating at 120 inches-per-second recorded all data channels. 
On-site verification of data was accomplished with a Honeywell 
1508b Visicorder. Post- processing of the data was through 
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TABLE II 



Instrumentation Specifications 



EQUIPMENT 
Strain gages 
Pressure gages 



TYPE 

CEA-350 ohms 
.25'' Tourmaline 



MAX RANGE 

50,000 microstrains 

10 kpsi, 97% Response 
Ratio 



Accelerometers PCB Piezoelectric 2500 and 100,000 G's 
Amplifiers Ektron 563FJ 



an HP-5451C Fourier Analyzer located at NPS . A schematic of 
the instrumentation is shown in Figure 4.1. To assist in 
data reviewing, an ingenious triggering mechanism was devised 
by electronic technician Tom Christian. This triggering 
device consisted of a small DC circuit wrapped around the 
charge. At time of detonation the wire breaks registering 
a zero voltage reading on the recorder. 

Instrumentation wires ran down from the connection box 
on the test structure, along the bay's floor to the WCSF ' s 
instrumentation shack. Figure 4.2 is a photo of the Honey- 
well recorders and their associated amplifiers in the 
electronic's shack. It should be noted that a redesigned 
connection box was constructed and implemented for this shot. 
On the previous UNDEX experiments a brass electrical box was 
installed, which did not afford adequate protection to the 
connections. To remedy this weakness, a new box was fabri- 
cated out of 1.25 inch, A36 structural steel plate. 
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Figure 4.1 Instrumentation Schematic 
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Figure 4.2 Honeywell Recorder and Amplifiers 




B . SHOT DAY 



Thursday, 6 December 1984, was set as the shot date. A 

■rtf 

drawing of the test rig is provided by Figure 4.3. The 




Figure 4.3 Test Rig Model 

weather was ideal for an UNDEX experiment--clear and sunny 
skies, calm seas and no wind. After re-examining the equip- 
ment and performing calibration checks on the sensors, the 
riggers took over preparing the test structure for the bay. 
Figures 4.4, 4.5, and 4.6 show the shot sequence. Immediately 
after the explosion, everyone sensed something had gone wrong 
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Figure 4.4 Crane Hoisting Test Panel 
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Figure 4.5 Test Panel Submerged/Boom Removal 
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Figure 4.6 First Explosion 
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with the shot as the water plume was smaller than expected, 
and the explosive sound was muffled. Tom Christian, stationed 
in the instrumentation shack, yelled out moments later that 
he had picked up some very weak signals on all channels, but 
was unsure of their meaning. Figure 4.7 reveals the test 
structure as it is being hauled out of the water following 
the explosion. No significant external damage was apparent 
and it was noticed that the free field pressure transducer 
was intact and still properly oriented. It was not until the 
test structure was repositioned on dock, that a large amount 
of unexploded TNT was discovered on top of the plate, see 
Figure 4.8. The appearance of the TNT on the plate was con- 
clusive that the full eight-pounds of explosive did not go 
off. The unexploded residue was collected and weighed approxi- 
mately one and a half pounds. During post-analysis, by 
examining the pressure graphs from the Fourier Analyzer, it 
was determined that approximately 0.16 pounds of TNT actually 
exploded. This was obtained by inputting the recorded value 
of 820 psi into equation 2.3 and solving for charge weight. 

The plate surface was carefully inspected with no visible 
damages noted. Deflection measurements were not taken as a 
micrometer was not available. The electrical connection box 
remained water-tight and all but one of the sensors were 
still operational. After a short deliberation and determin- 
ing the availability of another eight-pound charge, it was 
decided to fire another shot later that same day. WCSF 
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Figure 4 . 7 



Test Structure After First 



Explos ion 
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Figure 4.8 Unexploded TNT from the First Shot 





personnel had stated that this was the first incomplete 
explosion of this type they had experienced. Incidentally, 
it should be pointed out that firing a small primer charge 
before a larger explosion is not an uncommon practice. 
Experience has shown that a primer tends to aid in gage 
calibration and loosens up residual plate stresses. 

Approximately three hours later, another explosion 
occurred. This time two blasting caps were used to ensure 
full detonation of the TNT. It was apparent from the noise 
and water plume that the second shot was significantly 
larger than the early one. It was later determined from 
pressure data graphs that the peak pressure pulse experienced 
by the plate for this second shot was 4217 psi, which is 
equivalent to a 9.84 pound TNT charge. Figure 4.9 shows the 
water plume produced by the second shot. One of the pneumatic 
bladders ruptured from the force of the explosion, as shown 
in Figure 4.10. Another post shot view. Figure 4.11, shows 
the damage experienced by the test rig. Notice how two of 
the test platform support cables have been torn free. 

A view of the underside of the test structure and cable 
runs going into the electrical connection box is shown in 
Figure 4.12. In the past the electrical cables have been 
pulled loose from the box, exposing the connections to the 
environment. Following this second explosion, Tom Christian 
joyfully announced that data had been collected. 



44 




Figure 4.9 



Water Plume from the Second Shot 
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Figure 4.10 



Ruptured Pneumatic Bladder 
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Figure 4.12 



Clamped Cable Run into Connection Box 
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V. RESULTS AND RECOMMENDATIONS 



A. PHYSICAL RESULTS 

As mentioned in Chapter IV, the first shot did not pro- 
duce any noticeable physical results. However, the second 
shot did achieve large deformations as evidenced by Figures 
5.1 and 5.2. Both stiffener's ends sheared and there are 
tears along three of the sides. The exact tear lengths are 
recorded on Figure 5.3. Not visible in any of the photos is 
the fact that several strain gages became detached from the 
plate surface during the testing. Additionally, a few of the 
strain gage lead wires were severed. This severing was 
probably due to a guide wire shackle that parted from the 
ruptured flotation device following the shot. 

Dial indicator deflection readings were taken while the 
aluminum plate was still clamped to the test structure, and 
are presented in Figure 5.3. An EPSA run was made using a 
9.84 pound TNT charge to determine predicted deflections. 

EPSA predicted and actual displacements are compared on 
Figure 5.3. There are large discrepancies noticed in these 
comparisons, especially near locations where plate tears 
exist. EPSA does not take into account possible secondary 
effects which may have loaded the plate, most notably the 
loading generated by water particles striking the plate due 
to the downward rigid body action of the test structure. 
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Figure 5.1 Shearing at Both Stiffener Ends 
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Figure 5.2 Large Plate Tear at the Side 
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Figure 5.3 Displacement-EPSA 9.84 lbs TNT Predicted vs Actual 



EPSA assumes a perfectly symmetric loading pattern, without 
any stress concentrators and does not incorporate fracture 
mechanics . 

B. DATA INTERPRETATION 

The majority of data interpretation was conducted at NPS 
using the HP-5451C Fourier Analyzer. Sensor data recorded 
on FM tape was transferred from the Honeywell recorder to a 
5451C disk. All channel data was recorded at 120 inches-per 
second ( ips ) and played back to the disk at 1.875 ips . A 
real time conversion factor of 64 is necessary to change from 
120 to 1.875 ips. Time zero for data reduction purposes is 
the exact moment of detonation. The triggering channel was 
instrumental and served as a reference in transferring data 
from the recorder to the disk. All seventeen channels are 
permanently stored both on tape and disk at NPS. Each data 
channel consists of fifteen storage records, with each record 
tracking 20.47 msec at 1.875 ips time or 0.319 msec real time. 

Table III gives a summary of data, storage locations and 
conversion scale factors used for shots one and two. Sample 
digitized output data graphs are provided in Appendix A. 

Time was not available to write a computer program to change 
the scales on the axes to reflect proper units or actual 
times. However, reference values are indicated on each 
graph . 

As mentioned earlier, the digitized incident pressure 
data from shot one yielded a peak pressure of 820 psi, 
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TABLE III (CONTINUED) 





u 


o 


u 


u 


o 




o 


o 












M 


CD 


CL) 


CD 


d) 


a) 




CD 


a) 












Eh 


CO 


c o 


CO 


CO 


co 




CO 


co 














a 


a 


a 


a 


a 




E 


E 












i-Q 




























g; 


co 


i — i 


LO 


1 — 1 


o 




cn 














> 














co 


co 












H 














• 


• 












Dh 


CM 


CM 


CM 


CM 


CM 




CM 


CM 


CD 




























a 










< 


















cd 




























MH 




























U 




























0 




PQ 
























CO 
































U 










co 


co 


CO 


CO 




CO 


CO 


CO 




0 










G 


g 


G 


G 




G 


G 


- 














•H 


•H 


•H 


•H 




•H 


•H 


0) 




< 










Cd 


cd 


cd 


cd 




cd 


cd 


-p 














P 


u 


U 


u 




u 


u 


cd 




u 






CO 


O 




-p 


+3 


-p 




■p 


-p 


i — i 




CD 






W 


2 


CO 


CO 


CO 


co 




CO 


CO 


Ph 




n 






D 


M 


01 




ZL 


G- 




H 


11 






H 








Q 






JX 










G 




0 






<: 


< 






CO 










0 




o 






> 


PQ 




CM 








r- 


LO 






CD 








PP 




r- 




CO 




co 




CD 




pp; 




CO 








• 


, — | 


CM 




• 


• 


0^ 








-H 


< 


Q 
















cd 




i — i 






W 


< 


^ N 


.» s 


x— *S 


X-N 






X-N 


a 




i — i 






PH 


PP 


> 


> 


> 


> 




> 


> 






CD 




o 




















CD 








1 — 1 






r- 


CO 




LO 




r- 


LO 


U 




>i 




LO 






pH 


cn 


G\ 


00 




co 




0 




CD 










t 






CM 




o 


o 


CO 




G 




LO 






i — 1 


• 


• 


• 




• 


• 


CO 




O 




PH 








— 


— 






— ' 


' — - 


CD 




PC 




X 




















U 




























PL, 




II 




G 




























0 




















i — 1 




=8= 
























cd 


CD 






CO 




















■P 


cp 


PQ 




G 




















0 


cd 






0 




















Eh 


a 


6 




*H 


w 


























-P 


PI 


















II 


CD 






cd 


< 




















u 






u 


u 


















, — > 


0 


U 




0 


CO 


















-P 


CO 


CD 




r — 1 




















— ' 


CO 


o 
























PH 


CD 


a 




CD 


o 












a 








u 


0 




CP 


M 


o 


o 


LO 


r- 


o 


pq 


o 


o 




PH 






cd 


CO 








00 


o 


Eh 


1 — 1 


1 — 1 


U 








u 








• 


• 


1 — 1 


U 






CD 


T3 


CD 




O 


w 






i— i 


CM 




PQ 






-P 


i — 1 


CT^ 




-P 


> 












P3 






CD 


CD 


cd 




CO 














hQ 






a 


•H 


a 






o 












O 






0 


Ph 






cd 


o 












o 






P 




G 




-P 




















CD 


CD 


*H 




cd 














< 






i — 1 


CD 


cd 




HG 














Eh 






CD 


H 


u 








o 


LD 


LO 


O 


LO 


< 


o 


o 


O 


Ph 


V 


-p 


II 


Q 


o 


1 — 1 


r- 


<T\ 


o 


Q 


CM 


LO 


O 




CO 


i — i 




w 


ro 


CO 


co 


CO 






H 4 




< 


II 




0 


Q 


os; 


1 


1 


1 


| 


l 


O 


1 


I 






II 


> 


PQ 


o 


CO 


1— 1 


r— 1 


CO 


• — i 




CO 


CO 


II 








X 


Eh 


00 


o 


CO 


r- 


o\ 




o 


CO 




MH 


=8= 


II 


o 


CO 


CM 


CO 


CO 


co 


co 








- 


' — 


a 




Eh 






















PH 


CO 


> 


CO 




CO 


r- 






o 


rH 


CM 


co 












PP 


U 


O 


00 


<y\ 


1 — 1 


r— 1 


rH 


i — i 


• • 










O 


CO 


CO 


o 


o 


a 


a 


a 


a 


CO 










CO 






CO 


CO 


CO 


CO 


CO 


CO 


a) 




• 






2 


1 — 1 


CM 














-p 










PQ 


r— 1 


rH 


CM 


CO 




LO 


CO 


r- 


0 










CO 


< 


< 


PQ 


PQ 


PQ 


PQ 


PQ 


PQ 













55 



equating to 0.16 pounds of TNT at a real time of 2.05 msec. 
This value was compared to that of the plate's surface total 
pressure gage of 585 psi at time 2.13 msec. It is difficult 
to say why the peak surface pressure is less than the peak 
incident pressure but perhaps the time difference of 0.079 
msec accounts for this degree of decay. This statement is 
predicated on the assumption that both pressure transducers 
were the same distance away from the charge. Selecting the 
value of 0.16 pounds of TNT, another EPSA computer run was 
made to determine the impact of this weight. Examining the 
EPSA output indicates the possibility that some permanent 
plate deflections may have resulted from the first shot. 
Maximum center plate deflection was calculated to be 0.0147 
inches. Using the aluminum plate thickness of 0.25 inches, 
a ratio of deflection to plate thickness of 4% is attained. 
This is considered a small deflection and serves as a basis 
to ignore the impact of the first explosion. Another test 
performed to satisfy the hypothesis that the first shot did 
not have an overpowering effect on the test plate was a 
comparison of peak strains to yield strain. The typical 
yield strain for 6061-T6 aluminum is 4000 microstrains. The 
EPSA strains produced from a 0.016 pound TNT load are all 
well below yield. This indicates that the plate did not 
enter its plastic region. 

Peak strain comparisons are provided in Figure 5.4. The 
values within the parentheses, are EPSA predicted and the 
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4 7k (37. 5k) 
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Figure 5.4 Comparisons of Predicted to Actual Strains 



•ones outside represent experimental. It was disappointing 
to find such large discrepancies between gages. The measured 
peak values are haphazard, without any regard for consistency 
in reading higher or lower than expected. Two of the strain 
gages gave no output due to severed lead wires, and one indi- 
cated a reading beyond gage limits. Accelerometer magnitudes 
were equally useless. Both transducers registered values 
exceeding their maximum range. Sensor symmetry responses were 
investigated by comparing the outputs from three sets of 
strain gages and both accelerometers. Table IV gives the 
highlights of these comparisons. As alluded to earlier, the 
peak strain values are inconsistent but loading times are 
very close. Poor peak value comparisons add a great deal of 
uncertainty about the validity of the strain gage measurements. 



TABLE IV 



Symmetry Sensor Comparison 



SENSOR 



SHOT #1 



SHOT #2 



SG1&7 



No comparisons 

due to bad SGI values 



Arrival Times within 0.7% 
Peaks off by factor of 2.2 



SG3&6 



Arrival times within 1.0% Arrival times within 0.7% 
Peaks off by factor 2.1 no comparison, bad SG6 value 



SG10&13 not recorded 



Arrival times within 2.4% 
Peaks off by factor of 62 



X" 



Arrival times within 0.3% Arrival times within 2.7% 
Peaks off by factor 1.16 Peaks off by factor 1.4 



Notes: Peaks represent strains for the strain gages 

and G's for the Accelerometers (X") 



58 



The data pressure histories from -shot two, enclosed in 
Appendix A, seem to be the most useful data obtained from 
this experiment. The peak surface total pressure sensor 
value of 4.922 kpsi is larger than the free field peak of 
4.217 kpsi. Theoretically, plate total pressure should be 
twice as large as free field; but this is assuming an ideal 
rigid body response. An aluminum plate does not reflect 
perfectly, so any magnification between one and two is 
practical . 

There does exist a time difference between pressure 
transducers. Prior to the structure's water entry, there was 
exactly a nine foot charge to transducer separation distance. 
A time comparison attempt was made to check whether these 
distances remained the same once submerged. The results 
of this comparison are shown in Table V. Theoretically, the 



TABLE V 

Pressure Sensor Time Comparisons 

SHOT #1 SHOT #2 



SENSOR 


ARRIVAL 


PEAK 




ARRIVAL 


PEAK 




P(t) 
P(f ) 


2.116 

1.950 


msec 

msec 


2.132 

2.024 


msec 

msec 


2.390 

2.185 


msec 

msec 


2.403 

2.213 


msec 

msec 


Time change 
Dis tance 


0 . 166 
0 . 83 


msec 

ft 


0.109 

0.54 


msec 

ft 


0.203 

1.02 


msec 

ft 


0 . 190 
0.95 


msec 

ft 



Note: Distance is computed assuming the speed of sound, 

5 ft/msec 

shock wave should have reached the plate at 1.8 msec, as 
computed from equation 2.2. The free field pressure 
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transducer, although identical to the other, does consistently 
give earlier signal readings. Table V indicates that a 
stand-off distance of exactly nine feet may not have been 
maintained between the charge and pressure transducers. 

C. POST-SHOT EPSA TESTING 

Initial comparison of EPSA predictions and experimental 
results indicated that the EPSA model was inaccurate or had 
some major prediction weaknesses. However, there are too 
many variables present in this UNDEX experiment to use such 
a gross generalization. The mechanics incorporated in EPSA 
are sound and based upon singularly proven theories. During 
an experiment, although ideal conditions are sought, they 
are seldom totally attained. Adopting a more realistic 
approach, further tests were conducted using EPSA. The main 
goal of the post-EPSA tests were to try to achieve the measured 
deflections using the code. Computer printouts of data re- 
sults and PATRAN-G color slides were taken of the various 
EPSA runs , but they are too lengthy to include in this re- 
port . Appendix B provides sample EPSA input decks used for 
the different runs, should there exist an interest in 
recreating the data. 

The following describes a list of post-EPSA tests that 
were conducted to try to imitate test results. First a 9.64 
pound TNT run was made but the deflections were not signifi- 
cantly different from that of an eight pounder. Accounting 
for the fact that plate tears did occur, fundamental fracture 
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mechanics were attempted to be incorporated into EPSA. A 
change to the subroutine SM STIF allowed for modelling the 
plate with the stiffener ends cut. This provided a less 
rigid test panel but deflections remained symmetric. The 
stiffener ends were reattached and boundary conditions were 
written into the EPSA input deck to simulate actual plate 
tears along the top and side edges. The results from these 
cuts were quite interesting as this was the first time non- 
symmetric EPSA deflection patterns were produced. However, 
these cuts did not portray the time history or sequence 
whereby these tears developed. Finally, using data points 
from the second shot's classic free field pressure plot, 
actual pressure-time history information was inserted into a 
data deck. The pressure-time history run ignored charge 
weight and tearing conditions. The results were similar to 
a 9.84 pound charge. This was expected as it was from this 
graph that this weight was calculated. Using Table VI and 
Figure 5.5, a complete picture of EPSA deflection results 
are given. The term 'cut' in Table VI refers to the tear 
conditions observed from the second shot. None of the above 
steps reproduced the experimental results, but they did 
demonstrate some of EPSA's capabilities. 

D. CONCLUSIONS AND RECOMMENDATIONS 

After reviewing the results of this experiment and the 
others in this series, it is clear that stiffener tripping 
is not occurring. Stiffener tripping is not experienced due 
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TABLE VI 



Comparison of Actual 


Deflect 


ion to 


EPSA Predictions 














CUT 


CUT 


CUT 


^RK 


ACTUAL 


8# 


9.84# 


0.16# 


P-T 


8# 


9.84# 


P-T 


1 


.212 


.061 


.072 


.004 


.078 


. 143 


.163 


. 700 


2 


. 583 


.328 


. 366 


.031 


. 406 


.618 


.681 


. 679 


3 


1.195 


. 619 


. 676 


.063 


. 688 


.958 


1.030 


1.030 


4 


. 757 


. 328 


.366 


.031 


.406 


- 


- 


- 


5 


. 211 


.061 


.072 


.004 


.078 


- 


- 


- 


6 


.448 


.138 


.156 


.017 


.115 


.294 


. 328 


.329 


7 


.939 


. 311 


.341 


.047 


.317 


. 60 8 


. 661 


.657 


8 


1.368 


.353 


.385 


.050 


.396 


1.150 


1.240 


1.220 


9 


1.042 


.311 


. 341 


.047 


.317 


- 


- 


- 


10 


.524 


. 138 


. 156 


.017 


.115 


- 


- 


- 


11 


. 514 


.072 


.085 


.005 


.085 


.480 


. 519 


. 50 5 


12 


1.008 


.323 


. 361 


.031 


.395 


. 642 


. 707 


.705 


13 


1.490 


. 601 


. 654 


.063 


.667 


1.010 


1.090 


1.090 


14 


1.082 


.323 


. 361 


.031 


. 395 


- 


- 


- 


15 


. 514 


.072 


.085 


.005 


.085 


- 


- 


- 


16 


.474 


.061 


.072 


.004 


.078 


.438 


.484 


.473 


17 


.993 


.328 


.366 


.031 


.406 


.618 


.681 


. 679 


18 


1.424 


.619 


.676 


.063 


. 688 


.958 


1.030 


1.030 


19 


1.025 


.328 


.366 


.031 


. 406 


- 


- 


- 


20 


.590 


.061 


.072 


.004 


.078 


- 


- 


- 


21 


.314 


.138 


.156 


.017 


.115 


.133 


.152 


.159 


22 


. 664 


.311 


.341 


.047 


.317 


.360 


.400 


.408 


23 


.944 


. 353 


.385 


.050 


.396 


.515 


. 556 


. 551 


24 


. 713 


.311 


.341 


.047 


. 317 


- 


- 


- 


25 


.383 


.138 


.156 


.017 


.115 


- 


- 


- 


26 


.106 


.099 


.114 


.011 


.088 


.097 


.112 


.122 


27 


. 241 


.174 


.194 


.024 


. 187 


.172 


. 195 


. 206 


28 


. 334 


.172 


.194 


.024 


.194 


. 236 


.256 


.258 


29 


.245 


.174 


.194 


.024 


.187 


- 


- 


- 


30 


.129 


.099 


.114 


.011 


.088 


- 


- 


- 


NOTE: All 


deflections 


have units of 


inches 








The 


mark number 


refers 


to the 


locations in 


Figure 5 
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Figure 5.5 Deflection Locations 



to the shearing of both ends of the stiffener. Numerical 
analysis does not predict the shearing or tearing along the 
sides. From the data, the sequence of fractures is not 
discernible. Knowing the order of fracture occurrence 
would give a better clue to plate responses. Additionally, 
EPSA predictions do not accurately forecast experimental 
results. In general the code predictions are lower than 
measured . 

To observe stiffener tripping, the stiffener needs to 
have a compression mode throughout its entire length. 

Curvature must be introduced to the test panel. It is 
recommended for future studies to perform UNDEX testing 
using a stiffened shell model. This would provide all- 
around compression on the stiffener, promoting stiffener 
tripping instead of shearing. Also, it is necessary to 
incorporate tearing effects or fracture mechanics into the 
numerical analysis scheme. Symmetry loading conditions are 
not always observed from an UNDEX experiment and these 
need to be taken into account. 

It is recommended that UNDEX testing be continued at NPS 
in order to improve experimental techniques and to acquire 
stiffener tripping data. Once tripping has been accomplished, 
then investigations can branch out to the gross shell 
responses of T and Z stiffeners. Improved stiffener tech- 
nology would lead to a significant advancement in ship 
design . 
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APPENDIX A 



HP-54 51C FOURIER ANALYZER DATA HISTORIES 
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Figure A.l Free Field Pressure Gage--Shot #1 



10 " 
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Figure A. 2 Surface Total Pressure Gage — Shot #1 



10 " 
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Figure A. 3 Free Field Pressure Gage--Shot #2 



o 




O l I I I 
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Figure A. 4 Surface Total Pressure Gage — Shot #2 



appendix b 



EPSA SAMPLE DATA DECKS 

JL 0 L 4 T F 

2500 1 0 0 0 *5 0 0 10 10 
. 00000 1 1 ? 5 0 0 25 

1 1 1 U 15 15 

1 5 25 c 1 1 0 4 24 40 50 1 1 0 

0 0 

10000000. .5 .000 25 5 40 0 00. . 25 0 0. . uOOQl 

’ S T I F • 1 1 u 

- 1 1 2 
4 16 0 

1 0 0 0. 0. 

1 6 1 o 4. 0 . 

-16 16 25 0 .45 

0 / 

575 1 1 ’5 25 

111611110 

2 1 ?6 1 0 0 1 0 

511611110 
41261 1 1 10 

11110.0.0. 

0 . / 

161 3 1 6 26 5 16 15 2 16 15 3/ 

1111121141151511 
15121514151551515132 
3134313581 1*12814 

8151711 72174175 

331 332334335/ 

151 1081 10157106710151210 

15 l q 10 1525 10 1 70872 1570 
11311151311 I5l'44l5 18l5l q 

81 881 315781579878 

87° 15 12 s 15 12 9 15 108 15 l 9 q 

15 25 8 15 25 9 177 17 109 77 

*71015771571001107110 

1 8 Q 88 Q 1588 )5 p 0 9 1 I 25 10/ 

1 6 23 3 a 2 5 5 3 *3 3 1 o 20 3 8 20 3 

320316)736173317316103 
810331031o7387 5373 

1643643343161331133 
923141231 4 133 ! 4 14315123 
151331514316123161437133/ 

.48 .260 1.250 

13 1 / 

0 / 

5.0 0.0000945 oOOOO.O 

1 106.166540 0.0554985 o.O 

8.0 29505.0 1.18 0.05* -0.1 a S 0.0012/ 

0/ 



Figure B.l 8 Pound EPSA Deck 
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JL a 1 E 

250 0 1 0 0 0 | 0.0 0 1 l 

. 0 0 0 0 0 1 / 

1 5 PS a 1 1 0 5 2 5 1 5 1 I o 
0 0 

10000 )00 . .5 . 00025" 4000 0 . .250 o. . ) o r. o 1 

' 3 T I F ' I 1 0 

-1 1 2 
4 I o 0 

1 0 0 a. o. 

lb 1 o 9 . o . 

-lb lb 25 0 .44 

0/ 

375 I 1 15^5 

1 1 I b 1 1 1 1 0 

2 1 2 b 1 0 0 i o 

3 1 1 b 1 1 t I 0 

4 1 2 b I l I I ' 1 

1 1 1 l 0 . 0 . () . 

0 . / 

1 b 1 3 1 b 1 2 1613216133112/ 

111112/ 

15 1 10 « 1 10 15 7 10/ 

3 4 3 6 4 3 1 b 4 3 3 7 3 H 73 

lb 7 3 3 13 33 133 1b 13 3 3 123 

3 12 3 | o 12 3 3 20 3 8 20 3 16 20 3/ 

.43 .250 1.250 

13 1 / 

0 / 

4.0 J.HM0 44 5 n 0 0,) 0.0 

2 104.166540 0.0554945 9.0 

7 

0.0 3046.4 >.00073 2343.6 0.0)3125 3230.44 0.io5b25 2342." 

0.02145 4 2 l *> . 0 4 0.02453 2^ 6o.8 0.0 3424 374H.8 
0.17548 250.0 

) / 



Figure B.2 Experimental Pressure — Time History 
Information 
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Jw PL4 I E 

2 5uO 1 > o 0 10 0 0 1 1 

.000001/ 

152571103111511) 

0 0 

1 0000000 . . 3 . 0 ,1025) 'loOOO. .1675 0 . 

• s r if - i i o 
- 1 12 
4 I 6 0 

1 0 0 o . u . 

15 1 0 0 . 0 . 

- 1 5 15 ? 5 0 . a-t 

0/ 

375 1 1 15 2 5 

111511110 

2 1 2 5 1 0 1 0 1 

3 15 1111 o 

3 5 1 5 0 0 ,i o | 

4 19 11110 

4 1 0 8 3 0 0 o 0 1 
42455111(0 
11110. I. 0. 

0 . / 

151110121513/ 

15 1 l / 

15 1 7/ 

3 4 3 6 431^^3373-^75 
1 5 7 3 3 1 3 3 « 1 3 3 l o 1 3 3 3 1 2 3 

5 1 ? 3 1 o 1 2 3 3 ? 0 3 * 20 3 1 r> 2i 3/ 



.46 .550 1.25 




13 1 / 




0/ 




9.0 ) . 0 u o u 5 a 5 


huoy.i . ) 


1 1 o 6 . 1 o o 5 4 1 


'» . 0 5 9 4 0 65 


9.8 4 ?P-, 1 1.0 


I.M " . 


0/ 





.ooo o i 



0 . O'M 2/ 



Figure B.3 9.84 Pounds, Cut--Top and Sides 
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